Conductance of d-wave superconductor/normal metal/d-wave superconductor 

junctions 
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We develop a theory of the conductance of superconductor/normal metal/superconductor junc- 
tions in the case where the superconducting order parameter has d-wave symmetry. At low temper- 
ature the conductance is proportional to the square root of the inelastic electron relaxation time in 
the bulk of the superconductor. As a result it turns out to be much larger than the conductance of 
the normal part of the junction. 

PACS numbers: 74.45.+C, 74.50.+r, 74.20.Rp 



At small voltages V, the current / through a super- 
conductor/normal metal/superconductor (SNS) junction 
can be written as 



I = J(4>) + G(4>)v. 



(i) 



Here cj> is the phase difference between the superconduc- 
tors, and J{4>) is a periodic function of <j) with period 
2tt. The time dependence of <f> is given by the Josephson 
relation: 



dtp 
It 



2eV 



(2) 



The first term in Eq. ([T]), describing the Josephson 
current, has been the subject of intensive experimental 
and theoretical studies over several decades, since the 
discovery of the Josephson effect. The second term, rep- 
resenting the dissipative current, has attracted relatively 
little attention, both on the theoretical and experimental 
sides. In the context of SNS junctions with s-wave sym- 
metry of the order parameter in the leads, this problem 
was considered theoretically in 0, [3, H, HI H, @, , and 
experimentally in Q . The common result of these works 
is that at low temperatures the conductance of the sys- 
tem is proportional to the inelastic relaxation time in the 
normal metal. A theory of conductance of a junction in 
the case when the superconducting order parameter has 
rf-wave symmetry has not been developed. In this paper 
we show that at low temperatures T the conductance, 
Gdnd, is proportional to the square root of the energy 
relaxation time in the bulk of superconductor, n n , and 
that it is much larger than the conductance Gn of the 
normal piece of the junction. 



The origin of the leading low-temperature contribution 
to Gdn d is similar to the Debye relaxation mechanism in 
dielectrics, or the Mandelstam-Leontovich mechanism for 
sound absorption in liquids with internal degrees of free- 
dom Due to the proximity effect the single particle 
density of states in the normal region, vn(s, 0), becomes 
e and (^-dependent. Here e is the energy of quasiparti- 
cle. According to Eq. ©, fjv(£, <fr(t)) changes in time. In 
the adiabatic approximation the electron population fol- 
lows the motion of the levels. As a result, the quasipar- 
ticle distribution function becomes nonequilibrium. Its 
relaxation leads to the entropy production, and therefore 
contributes to the conductance: 



TS = V 2 G 



DND ■ 



(3) 



The equation for the entropy production reads 



dr 



ds 



(l-/t 2 J 



Djj df df | {f-fthf 
v dxi dXj n n (r) 



(4) 

where f(e,r,t) is the quasiparticle distribution function, 
fth(s,T,t) — tanhe/2T is the equilibrium distribution 
function, n n is the inelastic relaxation time, Dij(e) is 
the diffusion coefficient, and v is the reduced density of 
states, measured in the units of the normal state density 
of states, vq. We assume that the latter is the same in 
both metals. 

The kinetic equation describing the dynamics of quasi- 
particles has the following form: 



V dt d Xi ij[T) dxi 



df 



v{f ~ hh) 
Tin(r) 



(5) 



The diffusion form of Eq. ([5]) is valid both inside the normal metal region if Ljv S> In, an d m the bulk of d- 



2 



wave superconductors at distances much larger than In- 
Here Ln is the length of the normal metal part of the 
junction, In,d are the the elastic mean free paths in the 
normal metal and d-wave superconductor, respectively. 

The last term in the left hand side of Eq. |(SJ) describes 
the generation of a nonequilibrium quasiparticle distri- 
bution function due to the motion of the energy levels. 
Qualitatively, the origin of this term can be illustrated 
using the example of a system of discrete Andreev levels 
labeled by their energy. Due to conservation of the num- 
ber of levels, their motion is described by a continuity 
equation 



dN(e) 

at 



de 



0, 



(6) 



where j e is the "current of levels" in energy space, and 
N(e) is the total density of states at energy e. If a level is 
occupied with a probability /, its motion causes a change 
in the occupation numbers of quasiparticles. The rate 
of change of the total number of quasiparticles at given 
energy, Nf, is determined by a continuity equation with 
a current fj s : 



d(Nf) , d(fj e ) 



dt 



de 



= 0. 



(7) 



Then using Eq. (jSJ) one can obtain Eq. |(SJ). 

In what follows, we assume that the last term in the 
left hand side of Eq. ©, describing the generation of a 
nonequilibrium distribution of quasiparticles, is nonzero 
only in the normal part of the junction. The reason is 
the following. At small T <C A , where A is the am- 
plitude of the order parameter, the voltage drop takes 
place mainly in the normal metal region. Indeed, the 
majority of quasiparticles incident on the normal metal- 
superconductor boundary experience Andreev reflection. 
A certain small fraction of the quasiparticles incident on 
the boundary in the direction parallel to the node in the 
quasiparticle spectrum in superconductor can penetrate 
the boundary. The part of the distribution function of 
these quasiparticles (the distribution function fx in the 
notation of Ref. [l(|), which is responsible for the im- 
balance of electron and hole populations, as well as for 
penetration of the electric field into the superconductor, 
decays away from the boundary on a length scale Le, 
which is of the order of the elastic mean free path in the 
superconductor, Le ~ Id <C Ln ll,[l2|]- This is the 
reason why we take into account only the electron-hole 
symmetric part of the distribution function (the distri- 
bution function / of Ref. 0) in Eqs. © and ©. 

Generally speaking, near the superconductor-normal 
metal boundary the system does not have d-wave symme- 
try. As a result, an s-component of the anomalous Green 
function is generated. At distances larger than ijv into 
the normal metal only this component survives and gives 
rise to the ^-dependence of the density of states. To de- 
scribe this effect in the normal metal piece of the junction 
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FIG. 1: Schematic picture of a d-wave superconductor 
(D)/normal metal (N)/d-wave superconductor junction. 



we use the Usadel equations for angles 9(e,r) and y(g, r) 
parameterizing the retarded Green function 13|, 14, 15 j 



le sin ft — 



D N d 2 9 



Dn 
4 



^| mi, 2^ = 0. 
or 



dr V dr 



(8) 



where is the diffusion coefficient in the normal metal, 
and we dropped the time derivative, which is justi- 
fied provided that eV -C E c = D^/L 2 N . The normal 
and anomalous retarded Green functions, and the lo- 
cal density of states in the metal can be expressed as 
g R (e,r) — cos0(e,r), f R (e,r) = — i sinft(e, r)e lx , and 
v(r,e) — Re cosft(e,r). Elastic electron scattering sup- 
presses cZ-wave superconductivity. Thus we consider the 
case £d <C Id, where £d is the superconductor zero tem- 
perature coherence length. 

Boundary conditions for Eqs. ((HJ) in the case of a con- 
tact between a <i-wave superconductor and a diffusive 
normal wire were derived in Ref. [161 ]. They depend on 
the angle a between the crystal axis of the superconduc- 
tor and the normal to the iVD-boundary (see Fig.[T]). In 
the case ai = an = 0, for a high transparency boundary, 
and at T,E C -C Ao, the boundary conditions reduce to 



6(e,x= ±L N /2) : 
X (e,x = ±L N /2) 



B = arcsnuV^-l) « 0.43, 



±- 



(9) 



In the case of s-wave superconductivity the system of 
Eqs. (O and ||5J) was derived in [nj. In the d-wave case, 
and for energies larger than the zero-energy scattering 
rate T (defined below), Eq. §5§ follows from the Boltz- 
mann kinetic equation for the quasiparticle distribution 
function. To see that Eq. is valid independently of 
the relation between e and T we illustrate its derivation 
in the Born limit of impurity scattering. To this end 
we use the quasiclassical equation for the Keldysh Green 
function g(e,t), which is a 4 x 4 matrix in the Keldysh 
and Gorkov-Nambu spaces [l(| , neglecting inelastic scat- 
tering: 

2^.5]+ + v F Vg - ^[er 3 + ir 2 A kl g} + ^-[(g),g]_ = Qp) 



3 



where [...,.. .] T denotes commutator and anticommuta- 
tor, T2 i 3 are the Pauli matrices in the Gorkov-Nambu 
space, Ak is a uniform d-wave order parameter corre- 
sponding to the direction k on the Fermi surface, td 
and vf are the normal state elastic mean free time and 



Fermi velocity in the d-wave superconductor. Introduc- 
ing the generalized distribution functions / anf f\ via 
9 K = (9 R - 9 A )f + (9 R T 3 - r 3 g A )h 0, Eq. (USD gives 
two equations for the latter: 



--((u)f-(uf)) + —((v)f-(vf)) 



U— + UVFVJl 

at 

11 HI 

uv F \7f - 2iA k wh = «u) /i - (u/i» + —({w) h ~ (wfi)), 



TD 



TD 



TD 



TD 



(11) 



where (. . .) stands for the average over directions on the 
Fermi surface, and the quantities u, v and w are defined 



-iTiT 2 (g H - g A )/4:, and 



by u = Tr(g R - g A )/4, 

w = — iTr T 2 (g R + g A ) / 4. The time derivative of f\ in the 
second of Eqs. (fTTj) can be neglected if we do not consider 
generation of branch imbalance in the superconductor. 
Expressing fx via vr^f using the second of Eqs. (fTTj) . 
and averaging the first one over direction of momentum, 
we obtain Eq. ([5]), with the following expressions for the 
diffusion coefficient and the density of states, valid in 
both Born and unitary limits: 



4ReiImi 



it 



t- 



2A k 



v D {e) 



Re- 



(12) 



(13) 



The quantity t(e) determines the retarded Green func- 
tion, g R = fat + i^Ak)/^ 2 
the following equation: 



A£, and is defined by 



t(s) 



t{e) 



2t d \y/W{e) 




(14) 



where (3 = ±1 corresponds to the Born and unitary lim- 
its of impurity scattering, respectively. The structure of 
Eqs. (fT2|) . (fT3|) . and (fl4|) is very similar to that appear- 
ing in the theory of the electronic thermal conductivity 
of (i-wave superconductors }17lll8|. 

Let us consider the geometry of the junction shown 
in Fig. [TJ We consider the ci-wave sup erconductor to 
be a stack of two-dimensional layers [19(. The planes of 
the layers are parallel to the axis of the junction. Wc 
also assume that within each layer the order parame- 
ter has d x 2_ y 2 symmetry, with order parameter given 
by Ak = Aocos[2(c£ — oll,,r)\ in the left and right 
leads, respectively. In this case the diffusion coefficient, 
Eq. dT3J) , for diffusion along the layers reduces to a scalar, 



Dij(e) = D D (s)S lj . Using Eqs. dHJ) and (JS} we get 



D D (e) 



^fV(e), N»r, 
el « r, 



(15) 



2ttA 



where 77(e) = 1 in the Born limit, and 77(e) « 
[2|e| ln(An/|e|)1/7rAn in the unitary limit. For the density 
of states we obtain [20( 



I ttAo 111 r 



el <r. 



(16) 



In the above equations the zero-energy scattering rate, 
r, is given by 



T = lmt(e 0) 



Aoe _7rAoTD , Born limit, 
75 — rr\ \j unitary limit. 

2t d ln(A T_D) ' •> 



(17) 

We now turn to Eq. (fTTj) . Formally, it should be sup- 
plemented with boundary conditions at the ND bound- 
aries. However, provided that the temperature is low 
enough, the time a quasiparticle spends in the normal 
region is small compared to T; m . Therefore, the third 
term in the left hand side of Eq. ([5]) can be substituted 
by lN(s,t)S(x), where 



In[s, t) — Ljy 



Of, 



lh 



ds 



,di/ N (e',t) 
dt ' 



(18) 



and ^/v(e, t) is the density of states in the normal part of 
the junction averaged over its length, 



1 

Ln 



L N /2 



L N /2 



dx vn{x, e, t). 



(19) 



Solution of Eqs. ([8]), (f9]), analogous to that presented 
in Refs. @, [5l| , yields the following e and f> dependence 
of the density of states in the normal metal region 

cos6» B , = 0, |e| < E c , 
;/ A = { sm0 B /0 B , (j) = vr, |e| < E c , (20) 
al, lei » B c . 
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The solution of Eq. |5} is 



/ - fth = lN7TR~ e 

ZVn 



\x\/L. r _ lD D T4, 
1 — 



(21) 



Substituting Eq. (J2TJ) into Eq. (4}, at eV < l/r irl we 
obtain the conductance 



Gdnd{.4>) ~ 2Gtv 



\EcTil 



Eq. (|22|) is the main result of this paper, which shows 
that the low-temperature conductance is proportional to 
the square root of the inelastic relaxation time in the bulk 
of the ci-wave superconductor. 

At low temperature the energy relaxation in e?-wave 
superconductors is determined by electron-phonon and 
electron-electron interactions. We would like to men- 
tion, however, that regardless of the relation between the 
electron-electron, r e _ e , and electron-phonon, r e _ p ^, re- 
laxation times, the conductance remains proportional to 
y/T e ~ph- The reason is that electron-electron scattering 
processes conserve energy. Thus if r e _ e <C T e _ p ^, at short 
time scales electron-electron scattering processes estab- 
lish an equilibrium form of the distribution function with 
a non-equilibrium value of temperature, which then re- 
laxes via electron-phonon processes. 

In the framework of BCS theory the electron-phonon 
relaxation rate in a d-wave superconductor at values of T 
that are not too small has the same order of magnitude 
as in the normal metal T e - P h oc O^/T 3 , where @d is the 
Debye energy. The reason is that the emission of phonons 
is not associated with a significant change in direction of 
the electron momentum. We note that in the context 
of high-T c materials the issue of the value and the T 
dependence of r in has not yet been settled [22| ■ 

According to Eqs. @ and (f2"2")l the conductance of the 
junction Gdnd (</>(£)) is a periodic function of cj> with pe- 
riod 2-7T, or a periodic function of time with period ir/eV. 
The amplitude of the oscillations is of the order of the 
average over the period of the oscillations conductance, 
Gdnd- The value of Gdnd depends on the ratio be- 
tween E c , T, and V: 

#s r>^ c >r, 



Gdnd ~ Gn&b V ^oTm < 



1 pE% 

»7(r) V t^t 




t, r > e c , 

; £ C >T>E 

£ c ,r>T. 



(23) 

Here the function 77(e) is defined after Eq. (fl~5j) . We also 
assumed that the normal state diffusion coefficients in 
the metals are the same. As a function of Q:l,r, Gdnd, 
Eq. (|22|) , has a maximum given by Eq. (|23p at ai = an = 



0, and vanishes at ol = an — 7r/4. In the latter case 
the conductance of the junction can be calculated in a 
way similar to calculations of the conductance of an ND 
junction, as in Ref. fl6| . and in that case Gdnd < Gjq- 
We note that Eqs. ^F2$ and (J23J) are valid at low 
voltages, eV < 1/T in . At eV > 1/T in , the same con- 
siderations lead to a non-analytic ^-dependence of the 
G D nd(V) cx V"- 1 / 2 , which can be obtained from Eq. (|23p 
by the substitution r in — > l/2eV. 

2 We would like to contrast our result, Eq. P3")) , with 
^) that for an SNS junction with s-wave pairing symmetry 
in the leads [J 0, 0, 0, !, 0, 0, H, in which the low- 
temperature conductance is proportional to the inelastic 
relaxation time in the normal metal, G$ns °^ T in . The 
origin of the difference with the s-wave case is that due 
to Andreev reflection the noncquilibrium quasiparticles 
cannot escape the normal region. 

The possibility of observing the effect discussed above 
is limited by the condition eV < l/rj„. At low temper- 
atures r in is long, and this condition is quite restrictive 
for junctions with leads made out of low-T c superconduc- 
tors. For junctions made out of high-T c superconductors 
the value of Tin can be much smaller, and the above re- 
quirement on voltages is much weaker. The values of the 
relaxation times in high-T c superconductors are not well 
established. If we estimate r\ n ~ lOps at T ~ 2QK in 
YBaCuO _23jj, we get the region of applicability of our 
theory V < lOO/xl/. We also note that the length of the 

junction is limited by the diffusion length DpfT^L e that 
is determined by the electron-electron relaxation time in 
the normal metal, T^L e . 

The obtained results are valid well below the critical 
temperature, T <C T c . Near T c , Andreev reflection at an 
ND boundary, which gives rise to the results presented 
above, becomes inefficient. Also the penetration length of 
the electric field into superconductor diverges as T — > T c , 
and the resistance of the junction can be calculated in the 
same way as a resistance of an NS boundary . 
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